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I We define correlational (von Neumann) entropy for an individual quan- 

t — , 

I turn state of a system whose time-independent hamiltonian contains random 



parameters and is treated as a member of a statistical ensemble. This entropy 
is representation independent and can be calculated as a trace functional of 



o 

I the density matrix which describes the system in its interaction with the noise 



source. We analyze perturbation theory in order to show the evolution from 



^ ' the pure state to the mixed one. Exactly solvable examples illustrate the 



use of correlational entropy as a measure of the degree of complexity in com- 
parison with other available suggestions such as basis-dependent information 
entropy. It is shown in particular that a harmonic oscillator in a uniform field 
of random strength comes to a quasithermal equilibrium; we discuss the re- 
lation between effective temperature and canonical equilibrium temperature. 
The notion of correlational entropy is applied to a realistic numerical cacula- 
tion in the framework of the nuclear shell model. In this system, which reveals 
generic signatures of quantum chaos, correlational entropy and information 
entropy calculated in the mean field basis display similar qualitative behavior. 



I. INTRODUCTION 



The concept of entropy is fundamental for many branches of physics and other sciences 
deahng with systems which reveal a certain degree of complexity and disorder. As stressed 
in the monograph , "entropy is not a single concept but rather a family of notions" . This 
monograph and the earlier review article contain historical information and give many 
examples of different approaches to the idea of entropy and numerous applications. 

In relation to quantum theory, the mainstream of development is formed by four main 
overlapping lines. They can be referred to as thermodynamical (Boltzmann - Gibbs) entropy, 
quantum ensemble (von Neumann) entropy, information (Shannon) entropy and dynamical 
(Kolmogorov - Sinai) entropy. Since the general description of a quantum system, including 
its interaction with the environment, time development and relaxation to equilibrium, can 
be given in terms of the density matrix the von Neumann definition seems to be the 
most fundamental. For a system in an equilibrium with a heat bath, the density matrix (and, 
accordingly, von Neumann entropy) is eqivalent to that in the canonical or grand canonical 
thermal ensemble. The evolution of a closed quantum many body (gas-like) system from 
a random initial state is shown to lead to the same values of macroscopic observables 
as for the thermal equilibrium described by the microcanonical ensemble which has a clear 
semiclassical limit as the equipopulation on the energy surface in phase space. The ensemble 
entropy cannot be represented as an expectation value of a dynamic variable expressed by 
an operator in Hilbert space. However, being a trace functional of the density matrix, it 
is invariant under unitary basis transformations. For a pure wave function, as that of a 
stationary state in an isolated system, basis-independent von Neumann entropy vanishes. It 
reflects the mixed character of the quantum state with incomplete information. 

Information entropy, with traditional applications in communication theory, is expressed 
in terms of probabilities rather than amplitudes. Therefore it is representation dependent 
being different for different choices of the set of mutually excluding events. In quantum 
systems, one can find information entropy of individual eigenstates with respect to a specific 
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basis. All correlations between the amplitudes of different components of the wave function 
are suppressed in this definition. Averaging information entropy over some ensemble of 
quantum states, one gets a measure of average complexity of those states. At this stage, 
the similarity between information entropy and von Neumann ensemble entropy can emerge 
if one can establish an appropriate correspondence between the ensembles used in the two 
approaches and the basis utilized in calculating information entropy . Thus, for canonical 
equilibrium thermal ensembles, the correlations are destroyed by the random interaction 
with the heat bath so that the density matrix is diagonal in the energy representation for 
the system. In this case, the eigenvalues of the density matrix give the occupancies of the 
stationary eigenstates of the isolated system which could be directly used for constructing 
information (= thermodynamical) entropy 

The so-called dynamical entropy, extensively studied during the last decade ||9Hl4 



IS 



essentially information entropy applied to a random sequence of measurements of quantum 
observables. Apart from the intrinsic complexity of the system under study, this construction 
reflects special features of the quantum measurement process. This entropy depends not only 
on the initial state but, in addition, on the observable and on the way of performing the 
measurement. It can be defined so that it give classical Kolmogorov - Sinai entropy in 
the corresponding limit of fine-grained phase space when the sequence of measurements may 



be described with the aid of symbolic dynamics |16 HT8|Jl^ . 

Information entropy used as a tool for quantifying the degree of complexity of individual 
quantum states |ll9|-|2^ shows delocalization of the wave function in a given basis. How- 



ever, as a rule one can find a basis, or a family of bases, which are singled out by physical 
arguments specific for each system. The delocalization length in such a representation man- 
ifests complex character of the state and can be quantitatively related to other signatures of 
quantum chaos. This exceptional role is naturally played by the coordinate representation in 
billiard-like cases and by the quasienergy basis in the problems with a periodic pertur- 
bation [^. For realistic many-body systems with strong interaction between constituents. 



the mean field represents the exceptional basis where the local correlations and fluctuations 



of adjacent stationary states are separated from their regular evolution along the spectrum 

m 



As shown in large scale nuclear shell model calculations p3|,[Y|], the representation de- 
pendence of information entropy might be considered in some respects as an advantage 
which provides a useful physical measure of mutual relationship between the eigenbasis of 
the hamiltonian and the representation basis. Moreover, chaotic dynamics make different 
states with close excitation energy and the same values of exact constants of motion "look 



the same" p4|, i.e. have similar observable properties. This is nothing but a microscopic 



picture of thermal equilibrium ||2^. After averaging over a narrow energy window in a high 



level density region, information entropy in the mean field basis becomes a smooth function 
of excitation energy and carries [^,0] the same thermodynamic contents as thermal entropy 
found for the microcanonical distribution from the level density. Being calculated in a ran- 
dom basis, the magnitude of information entropy of generic states in a complex system is 
typically on the level predicted by random matrix theory and does not display any regular 
evolution along the spectrum. 

The goal of the present paper is to explore the possibility of describing the degree of 
complexity of individual quantum states using the von Neumann definition of entropy and 
applying an external noise which converts a pure state into the mixed one. We do not con- 
sider the perturbation to be weak; therefore the resulting mixed state depends explicitly on 
the noise properties and gives a description of the ensemble "system plus noise" . There ex- 
ists a vast literature discussing the behavior of regular and chaotic systems under the change 
of parameters of the hamiltonian, see for example [p^^ -|37| . Multiple avoided crossings of the 



energy terms in a function of parameters reveal strong mixing and drive the system to the 
chaotic limit. The analogy of level dynamics with that of the one-dimensional gas of collid- 



ing particles is very productive for studying the spectral statistics p8| , |27|j28| , p^j39| , |i0[ ] . Here 
we assume that the ensemble of the parameter values is defined by a distribution function 
and calculate the density matrix and von Neumann entropy for a given energy term. Using 
exactly solvable models, we show essential features of represent at ion- independent entropy 



obtained according to this definition, its similarity to and distinction from information en- 
tropy. Even for the simplest systems as a harmonic oscillator in a random uniform field, the 
resulting steady states are far from trivial. We also give an example of a realistic numeri- 
cal calculation for a many-body system of fermions (a nucleus ^^Mg) which shows that our 
ensemble entropy ("correlational" entropy) is a smooth function of excitation energy and 
therefore may be used as a measure of the degree of complexity. 

II. DENSITY MATRIX AND CORRELATIONAL ENTROPY 

We consider a quantum system interacting with a surrounding. The interaction will be 
parameterized by a set of real parameters A in the hamiltonian, H = H{X). The energy 
spectrum of the system is assumed to be discrete. The eigenfunctions \a; \) of the system, as 
well as its energy levels Ea{X), evolve with A. For a complicated system, the level crossings 
are avoided so one can continuously follow these energy terms. 

At a fixed value of A, one can use any complete orthonormal basis \k) to study the 
evolution of the eigenstates in terms of the amplitudes C^(A), 

k 

Instead of the wave function (|l]), one can also use the density matrix p^"'^ whose elements 
are 

pS(A)=C,"(A)C,T(A). (2) 

p^"^ is a hermitian matrix in Hilbert space of the system. For a pure state |a;A), the 
descriptions in terms of the wave function ([l|) and the density matrix (Q) are fully equivalent. 
The obvious properties of the density matrix are the normalization 

Trp(")(A) = 1 (3) 

and the matrix identity 
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(p(")(A))' = p(")(A) (4) 

which shows that the eigenvalues of this matrix can be only or 1. Actually, the density 
matrix @ is diagonalized in the eigenbasis A). Only one eigenvalue, for the original state 
a' = a, is equal to 1 while the rest of the eigenstates belong to the degenerate subspace 
with the zero eigenvalues of p'^"-'(A) so that p*-°^(A) is the projection operator onto the state 
I a; A). Inversely, these properties can be taken as a signature of a density matrix describing 
a pure state. 

Now we assume that the interaction parameters A are random and have to be considered 
as members of an ensemble characterized by the normalized distribution function 'P(A), 
/ dW{\) = 1. Then the description in terms of a wave function becomes impossible, and 
our system is described by the density matrix (here and below the overline refers to ensemble 
averaged quantities) 

plt=WCW= JdXV{X)p'CM- (5) 

This is still a hermitian matrix with trace equal to 1. But, generally, the operator identity 
(|) is not valid anymore. The eigenvalues p^^^ of the matrix 

P^"V)=p1"V), (6) 
are nonnegative numbers between and 1, 

(Pl"^)'<pt"^- (7) 

These eigenvalues can be interpreted as mean occupation numbers of the eigenstates 
for a system which was brought into the contact with an external source being originally 
in the intrinsic state The exceptional case of a pure wave function is recovered for a 
fixed parameter, V{X) = 6{X — Xq). In notations of eq. and later on we distinguish the 
eigenstates \i>) of the density matrix from the eigenstates \a) of the hamiltonian by using the 
parentheses and the angular brackets, respectively; the dimensions of both sets are equal. 
The statistical distribution of occupancies can be characterized by von Neumann entropy 
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= -Tr { In p^") } = - E P^f^ In p[,"^ • (8) 

V 

This entropy, being still attributed to a single original energy term reflects correla- 
tional properties of the system subject to different levels of noise. Therefore we will call it 
correlational entropy although the definition is quite similar to that of standard ther- 
modynamic entropy in canonical thermal ensembles |]^. In contrast to information entropy 
of a given complicated state in a fixed basis |A;), which was used in the studies of quantum 
chaos [|19H21|J2|], 

r = -Y.'^t^^wt. w^ = {ci)\ (9) 

k 

correlational entropy (|^) is invariant and does not depend on the original basis of 
simple configurations. Obviously, this is a consequence of correlations between different 
components of the eigenfunction which are absent in the probabilistic definition of eq. (|^). 

For a pure quantum state |a;A), correlational entropy vanishes independently of the 
degree of complexity of the system. Thus, S*" characterizes the intrinsic term |a) as a 
member of the statistical ensemble. In general, S*" has an order of magnitude of InA^^Q, 
where is a number of the eigenstates \v) which have the occupancies p[f^ noticeably 
different from zero. The maximum possible value of correlational entropy is InA^ where iV 
is the dimension of Hilbert space. This value would correspond to the "microcanonical" 
density matrix, p[f^ = const = 1/N. Note that the information entropy has the order of 
magnitude of In where is a number of essential components in the stationary wave 
function |a; A) expressed in an original basis which was used in the definition of information 
entropy. Although formally the maximum value of information entropy is also InA^, its 
typical value in the Gaussian orthogonal ensemble of random real symmetric hamiltonians 
stays, due to the fluctuations, on the level of ln(0.48A^), see pT|j41| ,[7[| . In contrast to that, 



correlational entropy shows the degree of mixing, or decoherence, introduced by a given 
source of noise, regardless of the relationship between the resulting state and the original 
unperturbed basis. 
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In next two sections we consider simple examples which allow us to obtain exact solutions 
and shed some light on main properties and physical meaning of new entropy. 



III. PERTURBATION THEORY 

We start with the case of a narrow noise range AA that is small as compared to the scale 
of the parameter values which would lead to a considerable change of the wave functions. 
This is the perturbative regime. We can assume that the distribution function V{X) is 
concentrated near A = and find the perturbed wave function which starts its evolution for 
A 7^ from the unperturbed state |0). The state |0) acquires the admixtures of states \k ^ 0) 
which are given by standard perturbation theory (we assume the absence of degeneracy). 

With the hamiltonian H = Hq + XV where the perturbation V is an off-diagonal operator 
in the eigenbasis of Hq, the density matrix of the state |0) is, up to the second order in 
A, 

Plm = ^lO^mO I 1 ~ A^ ^ ^ I ~ A (SiQ—^ + 6mO — - 

V k^O J V em (-1 

I jy2 I ^O^mO I ^^o(^ ~ '^"'O) ^ K^fc^fcO I ^mo(l - 5iq) ^ VjfcVfcQ [ ^^^^ 

Here A and A^ are the mean values over an ensemble of noise, and the notation = Ek — Eo 
is used for the energy denominators. The density matrix ([10|) incorporates two effects, the 
redefinition of the original wave function |0) due to the perturbation XV, and the transition 
from the pure state to the density matrix. The first effect is the only one if the perturbation 
is fixed while the second effect appears because of the ensemble distribution of perturbations. 

It is easy to see that the first order correction to the density matrix does not change 
the purity condition (^, = p. The decoherence occurs only in the second order. The 
diagonalization problem for the matrix (0) can be solved to give, within a needed 
accuracy, two nonvanishing eigenvalues p^, as the roots of the characteristic equation 



p{l-p)-w{AXf = 0, (11) 
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and N — 2 zero eigenvalues. As expected, the statistical mixing is driven by the mean square 
fluctuation (AA)^ = — of the noise level. The decoherence rate is determined by the 
joint action of all virtual transitions, 

w = ^\vk\^, Vk = —. (12) 

k 

As seen from (0), only one state |0) has the eigenvalue po close to 1, 

Po = 1-«^(AAF- (13) 

The corresponding eigenfunction |0) = Z^fc has a large component ipj^^ ^ 1 — {w/2)X^ 

of the unperturbed state |0) and small admixtures of other unperturbed states, tp^^ ~ —Xvk, 
A; 7^ 0; the corrections are of higher order if the distribution function V{\) is even so that 
A = 0. 

In the approximation ([TlD, the new state |1) appears with a small but nonzero occupation 
factor 

pi = w{AXf. (14) 

The eigenfunction of this state is localized mainly in the subspace orthogonal to the unper- 
turbed state, ipit^ ~ VkW~^/'^, k ^ 0. The presence of noise removes the isotropic degeneracy 
of this (A^ — l)-dimensional subspace by singling out the direction of the multidimensional 
vector of transition amplitudes {vk}. Finally, the N — 2 degenerate states with the zero 
occupation factors pi, are orthogonal to the first two states. The high order perturbative 
corrections consequentially lift the remaining isotropy populating new combinations of orig- 
inal states. The decoherence process can be rather fast due to the added contributions of 
many distant admixtures so that perturbation theory can be valid at a very low noise level 
only. This "coherent decoherence" was discussed in a different context in [^. It is related 
to the selection of the most important (rainbow) diagrams in theory of disordered solids and 
in random matrix theory. 

The perturbed occupancies lead to nonzero entropy (||) 
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S = -polnpo- pilnpi^W{l-\nW), W = w{A\y. (15) 
The singularity at the origin imphes the infinite slope, dS/dW = — In TV, of growing entropy. 



IV. TWO-LEVEL SYSTEMS 
A. A two-level system with fluctuating level positions 

The simplest nontrivial but analytically solvable problem of a two-level system illustrates 
how correlational entropy is related to an interplay between the off-diagonal level mixing 
and their random diagonal displacements induced by the change of the parameter A within 
a range ~ A defined by an ensemble V{X). 

Here we consider two interacting states with a fiuctuating spacing so that the hamiltonian 
can be written, with the help of spin matrices, as 

H = ^{e-X)cr, + Va,. (16) 

The diagonalization of this hamiltonian at a fixed value of A is performed by the basis 
rotation through an angle ip/2 defined by 

2\/ e-A 
'^"^=A(A)' ^°^^=A(A)' 

where 



A(A) = ^(e-A)2 + 4V2 (18) 

is the level spacing. The amphtudes of the upper, (+), and of the lower, (— ), state can be 
written in the basis used in eq. (|T6|) as 



Cj = cos — , Ck = sm — , Ci = — sm — , Co = cos — . (19 
1 2 ^ 2 ^ 2 ^ 2 ^ ^ 

The instantaneous density matrix (0) for each of the two states, 

p(±) (A) = ^ [1 ± {cTx sin <f + a, cos , (20) 
10 



has the eigenvalues 1 and 0, and the corresponding entropy vanishes. After averaging over 
the ensemble of the values of A, we come to 

P^^^ = ^ [1 ± (^-^ + ^-c)] (21) 



where c = cosv? and s = siiup are averaged quantities which in generally do not satisfy 
+ = 1. This matrix has two nonzero eigenvalues, 

Pu = l{l + i^r), r = Vs^ + c2, u = ±l, (22) 

which are the same for both original states (±). Therefore both states possess the same 
correlational entropy 
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, 1 + T" 1 — T 

'\ + r) In h (1 — t) In 



(23) 



2 ' ' 2 

The exact result depends on the ensemble but, as a function of the averaging range A, the 
entropy value (|23|) evolves from 5 = at A = to S* = In 2 for the equipartition of the 
occupancies at large A when py = 1/2. 

As an example, we consider the behavior of correlational entropy for the ensemble with 
the uniform distribution of A in the interval (—A, A). With the definition (|l^) we have in 
this case 

e + A + A(-A) 1 A^ WAX! /^.x 

s = — In ^, c=— A-A - AA , 24 

A e-A + A(A) ' 2A ^ ^ ^ ^ ' ^ ^ 

The simplest case corresponds to the original degeneracy, e = 0, when 

r x/T+T^+l „ 2V 

s = - m , , c = 0, r = — — . (25 

2 yiTT^-i' 'A ^ ^ 

The averaged parameter s is small, s ~ rln(2/r), if the noise level is high and the intrinsic 
mixing is weak compared to the typical random level spacing, r <^ 1. At strong mixing, 
r ^ 1, s reaches 1. In this limit the levels are split by the dynamic interaction and the 
noise is ineffective in changing the population and reaching the decoherence. The entropy 
value correspondingly evolves from In 2 to 0. This behavior is opposite to the evolution of 
information entropy of the same states, 
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1 + COS 09 , 1 + COS 0? 1 — COS (Z3 , 1 — COS 09 

m m , (zDj 



2 2 2 2 

as a function of the mixing strength for a fixed diagonal spacing; J — > at weak mixing 
and J —> In 2 at strong mixing. In this example, information entropy measures merely the 
delocalization of the eigenstates with respect to the original basis. 

In Fig. 1 we show correlational entropy (^) as a function of r for different values of the 
ratio X = f/^- At x < 1 (Figs, la and lb correspond to x = 0.2 and 0.6, respectively), 



the results are qualitatively similar to those at e = 0, eqs. ( ]23|) and (pS]). With no dynamic 
mixing, r = 0, entropy is close to its maximum value, 5* ~ In 2—%^ at small x- As the relative 
strength r of the dynamic mixing increases, the value of entropy drops, 5* ~ [l + ln(4r)]/(4r) 
at large r. At x > 1, the entropy behavior changes since the noise level is not sufficient 
for covering the original level spacing, see Figs. Ic and Id, for x = 1 ^ind 5, respectively. 
In this situation, weak dynamic mixing increases entropy which starts from zero at r = 
and reaches its maximum value, much lower than In 2, at r = With the mixing strength 
increasing further, the level repulsion prevails which leads again to very low entropy. 

A complementary picture of Fig. 2 shows correlational entropy (|23|) as a function of the 
relative noise strength x^^ = A/e for four different values of the intrinsic dynamical mixing 
= 2V/e = t/x- Generally, entropy increases with x~^ and asymptotically approaches the 
maximum value In 2. For weak dynamical mixing (, the growth of correlational entropy is 
very steep in the vicinity of the point x = 1 where the accidental level crossing is possible 
for the first time, compare the change of behavior at this point in Fig. 1. The stronger is 
the dynamical mixing the smoother gets the increase of entropy. Entropy remains low as 
long as the intrinsic level repulsion governed by the parameter ( dominates. 



B. A two-level system with fluctuating mixing 

We consider the same hamiltonian (|1^) assuming now that the diagonal elements are 
fixed (we change the notation (e — A) — » e) but the off-diagonal coupling V fiuctuates. 
For definiteness, we assume that V is uniformly distributed between v > and —v. This 
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situation is simpler than that considered in the previous subsection because the only relevant 
parameter is that of perturbation theory, the ratio 

K = — 27 

e 

of the level spacing due to the mixing to the unperturbed spacing. 



The eigenvalues of the density matrix are still given by eq. (p2|) . Due to symmetry of 



perturbation, the average value of sin ip vanishes while the average value of cos ip is 

1 



c= - ln(K + V1 + k2). (28) 

With weak mixing, c ~ 1 — and entropy (|23|) is low, 

S ^ (kV12) ln(/tVl2), K < 1. (29) 

When the mixing intensity increases, c falls down, and entropy approaches its limiting value 
for equipopulated levels, 

5~ln2--^ln2-i(^) . ^ > i. (30) 

One can note that here the off-diagonal mixing plays a "natural" role of noise which creates 
entropy, in contrast to the previous example when it has been stabilizing the system by 
generating dynamic repulsion as a counterpoise to random noise; that situation can be 
considered as a prototype of the localization in disordered solids when the levels in the wells 
connected by tunneling fluctuate. 



C. Spin in a random field 

This two-state model is more general in the sense that both the mixing strength and the 
level spacing are fluctuating here. On the other hand, the level spacing in this model does 
not depend on noise so there is no level crossings. The hamiltonian describes a spin 1/2 
aligning along a random field. 
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H = a-n, (31) 

where the unit vector n has a random direction {9, if). Information entropy (^ of the 
eigenstates for a fixed field orientation, 

. If, l-cos^e , l + |cos^n 

does not depend on the azimuthal angle ^p. The density matrices (||) of the eigenstates |±) 
with the eigenenergies ±1 can be written as 

pW(n) = i(l±a.n) (33) 

and depend, contrary to (|3^), on both coordinate angles. 

After the averaging over the ensemble of angles, we find the eigenvalues of the density 
matrix (^) 

P, = \{1 + vt), v = ±l, (34) 



for any of the two original states (±). Here (compare eq. ([2^)) 



r = V = Jc^ + c=cos6', s± = sin6' exp(±iy9); < r < 1. (35) 



This allows one to find correlational entropy in the form (P3|) , 

5 = -i|lnl^+rlnl±^}. (36) 

The expressions ( p2D and (^) have identical structure. Correlational entropy reduces to 
information one in the case when the field direction is uniformly distributed on the surface 
of a cone with the fixed polar angle 9. The average density matrix is then diagonal in the 
original 2;-representation which is singled out by asimuthal symmetry of the noise distribu- 
tion. The two entropies coincide numerically at a point 9 = 9q where the effective angle 9o 
is defined by the condition cos^o = 
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D. Random matrix ensembles 



The 2x2 Hamiltonian matrix of an arbitrary two-level system can be presented as 

H=]^{u + a.w) (37) 

in terms of one scalar, u = TtH, and one vector, v = vn, parameters. The vector v is re- 
stricted to the (x, z) plane if time reversal invariance holds and becomes three-dimensional 
otherwise. Its magnitude v gives the level spacing. Both entropies / and S depend only on 
the unit vector n so that formulae (P^pB]) of the previous subsection remain valid indepen- 



dently of statistical properties of u and v. 

If all parameters of a hamiltonian matrix or some part of them can be random, the 
hamiltonian belongs to a random matrix ensemble. Gaussian orthogonal, GOE, and Gaus- 
sian unitary, GUE, ensembles of large random matrices attract special attention for modeling 
chaotic dynamics in quantum systems. In Gaussian ensembles, the hamiltonian matrix el- 
ements calculated in a fixed basis are taken as mutually independent Gaussian random 
variables. The idea of complete chaos implies that their distribution be invariant with re- 
spect to the choice of the representation basis \k). This invariance provides a rigid connection 
between variances of all matrix elements. However, it has been discovered that one gains ad- 
ditional possibility of describing special effects such as dynamical localization by permitting 
more flexible conditions for the variances. Some preferential basis there exists then where 
the localization takes place. This basis should be defined in each case in accordance with 
the concrete physical situation. Basis dependent information entropy or/and inverse par- 
ticipation ratio are used to describe the localization, see for example studies of the banded 



random matrix ensemble ||43|J4^ . Correlational entropy seems to present a universal and 
invariant characteristic of chaotic dynamics in presence of localization. We hope to return 
to this problem elsewhere. 

In Gaussian random matrix ensembles, all eigenstates have the same statistical distri- 
butions so that any of them can represent the generic features. Here we will only point 
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out the simplest properties of correlational entropy using the example of two-level Gaussian 
ensembles. According to (|35| - p6D , this entropy reaches its maximum value of In 2 under the 
condition of n = 0. This is obviously the case for invariant (orthogonal or unitary) ensembles 
where the vector n is distributed isotropically. For both the ensembles, mean information 
entropy, 

^ r difil{ip) = 2ln2-l, (38) 
Zn Jo 

or 

T^=l-Jdnii9) = ^, (39) 

is always lower than S. In eqs. (^8|) and (^9D the angular coordinates of the unit vector 
n are used, < (f < 271 and < < tt. The instantaneous quantity /(n) passes the 
maximum value of In 2 at the certain orientations of the vector n, as at the point 6 = tt/2 
in eq. (^2]), remaining smaller everywhere else. In fact, the relation / < S* is valid for the 
Gaussian ensembles of an arbitrary dimension N. As we have already mentioned, for ^ 1, 
correlational entropy is equal to S* = In A^ whereas mean information entropy reaches in this 
limit the value 7 = ln(0.48A^). 

At the same time, it can be shown by direct calculation that correlational entropy remains 
at the maximum value also when variances of all matrix elements have arbitrary values and 
the vector n ceases to be fully isotropic. Indeed, the probability distributions for ensembles 
of real symmetric matrices, 

cos^ + (T sm (f Zir 

and hermitian matrices, 

Vn^rMdn = ^ (41) 

q^^l cos2 61 + 0-2 Uf sin^ + ?| cos^ ipj sin^ 6 

remain even which yields n = as before. In eq. (^Tj), = {af + al) /2 whereas ai^2 and 
<^i^2 are variances of diagonal and off-diagonal matrix elements, respectively. 
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Correlational entropy drops to zero only when the off-diagonal matrix element vanishes 
identically. The randomness of interaction is of major importance. On the contrary, infor- 
mation entropy is sensitive to the ratios of variances which depend on the representation 
basis. 



V. HARMONIC OSCILLATOR IN A RANDOM FIELD 

Here we consider a problem with infinite Hilbert space, namely that of a linear harmonic 
oscillator in a uniform external field of a random strength. The diagonalization of the 
hamiltonian at any fixed parameter value is elementary but the diagonalization of the density 
matrix and calculation of correlational entropy are not trivial and we did not encounter the 



full solution in the literature although our results slightly overlap with those of ref. . 



A. The model 

At a given strength F, the hamiltonian of the problem, 

H = Ho-Fx = —p^ + -rmj^x -) 42 

describes the harmonic oscillator with the shifted equilibrium point, unchanged frequency 
and the energy spectrum 

En{\) = fiuj(n + 1/2 -\^) (43) 

where the integer n plays the role of the exact quantum numbers a in (|T]), and the dimen- 
sionless "noise" parameter 

F 

^ ^ {2hmuj^y/^ ^^^^ 

displaces the spectrum as a whole (no crossing in this model). 

The "natural" basis \k) of equidistant non-shifted states of a noiseless system, A = 0, is 
built with the annihilation and creation operators, 
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The shift of the equihbrium along the a;-axis to the point x = F/muj'^ = ^/2xo\ is the unitary 
t r ansf ormation 

D{X) = e-^v^^^op/^ = e^("^"'^) = D^-X) (46) 
with the obvious properties 

D{X)aD{-X) =a- X, D{X)a^ D{-X) = ~ X (47) 
so that (compare (^)) 

H{X) = D{X)HoD{-X) - X^huj. (48) 

The eigenvectors |n; A) at a given A are obtained from the original non-shifted states \n) = 
\n; A = 0) by the shift operator -D(A), 

\n;X) = D{X)\n). (49) 

Therefore the amplitudes of the eigenvectors (|l|) with respect to the noiseless basis \k) are 
given by the matrix elements 

C]:iX) = {k\D{X)\n) (50) 

of the shift operator p6|). The instantaneous eigenvector |n; A) carries an average number 
of the original quanta equal to 

n = {n; X\a^a\n; X) = n + A^. (51) 

B. Integral equation for the density matrix 

For the noise ensemble characterized by the (positively defined) distribution function 
V{X), the density matrix (|^) of the energy term |?t,) is 
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= / dXV{X){k\Di\)\n){n\D\\)\k') (52) 

where the unitarity of the D-transformation was used. This matrix represents the density 
operator 

pH = J dXVi\)D{\)\n){n\D\X) (53) 
in the noiseless basis. It is convenient to project instead the eigenvalue problem 

P^"V)=p1"V) (54) 

onto the set of states |n; A) with the noise A treated as the representation variable. Generally 
speaking, this set of states is not complete and maps the original problem on the subspace 
of eigenstates which belong to all nonvanishing eigenvalues p[i"^. In the new representation, 
the eigenf unctions 

cf>^-\X) = A^J^^.JviX){n\D\X)W), (55) 

where A^^^ are normalization constants, satisfy the integral equation 

/ dX'p'^''\X, X')<P^-\X') = p("Vl"^(A) (56) 

with the real symmetric kernel 

p(")(A, A') = ^/viX){n\Diy - X)\n)^v(y). (57) 

The diagonal elements of this kernel are simply the probabilities of the noise distribution, 

pH(A,A)=P(A). (58) 

They describe the ensemble itself and do not depend on the term n under consideration. 
The kernel (^) can be expressed as 

pW(A,A') = v'P(A)P(A')e-('-'')'/'L„,[(A-A')']. (59) 

in terms of the Laguerre polynomial defined by the series 
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, 79 I 

^0 [m\Y{n-m)\ 



It follows from (p3|-^) that the orthonormalized, (n;i/'|n;z/) = 6,^1^, subset of eigenvectors 



belonging to nonzero eigenvalues expressed in terms of the functions (^) as 

^) = (p1"0~'^' / dX,/v{X)<p'^:\X)\n; A) (61) 



if the mutually orthogonal solutions (^) of the integral equation ( ^6]) with the hermitian 
kernel (^7]) are normalized to unity which corresponds to the choice of A^J^^ = (^pjj*^) 

Starting with the completeness condition ^\ = 1; o^i^ easily comes to the 

kernel representation 

p(")(A,y) = Epi:"Vt"nA)0t"HA')- (62) 



In the original oscillator basis \k) = 0), the components of the eigenvectors (|55D are 

(eO, = (^1^;^) = [P^^V j d\^Wmk\D{\)\n). (63) 
Matrix elements of the shift operator are equal to 

{k\D{X)\n) = e-^^/^y||(-A)^-"Lr"(A^) (64) 

where 

is the associated Laguerre polynomial, L^{^) = Ln{^). Due to the symmetry property 

L'-iO = -i-O-'L^'iO, (66) 
n! 

the matrix ( |64|) is symmetric with respect to k and n. 

Up to this point, the results are valid for any energy term |n; A) and an arbitrary distri- 
bution function V{X). 
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C. Gaussian noise 



As an example of practical importance, we show here the exact solution for an oscillator 
originally in the ground state, n = 0, interacting with a source of Gaussian noise of a certain 
width A. The distribution function of the noise intensity is in this case 



and the ground state density matrix ( |^ reads in the A-representation (we omit the super- 
script n = 0) 

p(A,y) = ^e-(^^-^^'^)/^^^-(^-^')^/^ (68) 



Using the properties of the Hermite polynomials H,y{x), one can represent the kernel (pSQ as 

oo 

p{\,\') = Y.p,M^)M^') (69) 

in terms of the set of orthonormalized functions 

0.(A) = -^e-^^-^/^H.(v&;;A), (70) 

where 

= 2sinh(r//2)e-('^+^/2)^ (71) 
and the the parameter rj, < rj < oo, is defined by 



smhr]= — — — . (72) 

The representation (|69D obviously solves the eigenvalue problem (|56D for the ground state 
term. The normalization of the density matrix J2'^=oPi' = 1 can be checked directly. The 
value of correlational entropy defined by this density matrix is 

S = -Y,p,\np, = ^coth(r//2) -ln[2sinh(r7/2)]. (73) 



21 



The distribution function for the number of original quanta is given by the diagonal matrix 
elements of the density matrix in the unperturbed basis 



The average number of excited quanta equals 

k = Y.kfk = j dXX^nX) = F = A^, (75) 



k=0 

which leads to the simple value of average energy of the term 

E = 1i{pH) = huj{k + 1/2) -T^ = nu}{l/2 - A^), (76) 

in agreement with (^3]) and (|51|) for n = 0. The average energy of a given term is always 
lowered by the presence of noise; however, this energy cannot be attributed to the oscillator 
itself because it includes the contribution of the external noise source. 

According to (^) and (0), the mean number of original quanta excited by the external 
field can be written as 

^ (77) 



4sinh^(r7/2) 
Introducing A; as a new parameter, 



sinh I = l/{2\fk), coshr^ = 1 + l/(2fc), (78) 
we obtain the distribution of quanta (0) in the form 

When the external perturbation is turned off. A; — > and fk —>■ 5feo- In the opposite limit of 



strong noise, the mean number of excited quanta is large, k ^ 1, and eq. (|79|) reduces to 
the chi-square (Porter-Thomas) distribution. 



fk = {27rkk)-^/^ exp[-k/{2k)], A; > 1. (80) 
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We have to mention that, in the case of a fixed perturbation of a strength A, we would 
have a pure coherent state of a shifted oscillator with the Poisson distribution 

h = e-'^, (81) 

which reduces to the Gaussian distribution near the center k = k in the classical limit 
of ^ 1. The width of the coherent state in the occupation number representation is 
(A/c)^ = k and therefore the relative uncertainty of this number is small, Ak/k = k~^l'^. 
Contrary to that, in our random noise ensemble ( pOD the state is mixed and the fiuctuations 
of the number of quanta are not suppressed even in the classsical limit, 



[Aky = 2k\ (82) 

The situation is similar to that in the boson intensity interferometry (Hanbury Brown - 
Twiss) for many incoherent sources, see for example p6| . 



D. Relation to thermodynamics 

The eigenfunctions ( ffOD of the density matrix (|68|) in the A-representation formally coin- 
cide with the coordinate wave functions for a linear oscillator with the dimensionless 
"coordinate" ^ = v^smhr/A. The dimensionless hamiltonian for such an oscillator would be 

7^ = 1 (^2 _ ^2/^^2^ ^ (^tc, + 1/2) (83) 

where the annihilation and creation operators 

a = ^i^ + d/dO, a^ = ^i^-d/dO (84) 

are introduced. The energy spectrum of the operator (|83|), 6,^ = {u + 1/2), is equidistant 
and quantized in A-independent units at any noise magnitude. Since the density operator 
is diagonal in the basis of the eigenfunctions of (^3]), it is tempting to interpret our results 
as an equilibrium reached by the system under the infiuence of the noise. The equilibrated 
("dressed") system is represented by the effective oscillator (^). We can assign the physical 
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frequency uo to this oscillator and interpret the corresponding occupation number spectrum 
([n]) as the Planck formula with the effective temperature 



V 

Under this identification, the temperature scale is related to the range of noise, 

1 



(85) 



A 



2sinh{huj/2T)'' 
or to the mean square value of the random force. 



(86) 



(87) 



y/2 smh{huj/2T)' 

The density operator can be represented in terms of the effective oscillator defined above 



as 



p = exp i-hcuH/T) = p.W){u\, = 2 sinh l^—j exp [-jri'' + 1/2) j • (88) 

The effective oscillator is in the conventional thermodynamic equilibrium with a heat bath 
of temperature T, eq. (BSf). In an ordinary way, we can define its free energy 



J^ = Tln 



2 smh — 



(89) 



so that canonical thermodynamic entropy 



S 



df' 



(90) 



coincides with correlational entropy (|73D. Thermodynamic energy of the effective oscillator 
is equal to 



U = J- + TS = hu;{u + -) = —coth^ 

2 2 2T 



where the mean number of excited effective quanta is given by the Planck formula, 

1 



(91) 



exp{huj/T) — 1 



(92) 
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Let us consider an arbitrary dynamical variable of the original oscillator described by 
the operator 0{a^,a), 

= J dXV{X){X\0{a\a))\X) = J dXV{X)0{X, X), (93) 

where the operators and a of original quanta are supposed to be normally ordered with 
respect to the unperturbed vacuum; in the second equality ( p3l) the properties of the coherent 
state I A) were used. Since the parameters of the effective oscillator depend on temperature 
(^) via the factor i/sinh rj, the mean value of any operator (|93|) contains an additional 
nontrivial temperature dependence. Indeed, using the relations (|58| ) and (|69D , 



V{X) = p(A, A) = ^ p,0,(A)0,(A), (94) 



u=0 



we can express the mean value (^) in terms of the effective oscillator 



Tr 



pO 



u=0 

a'' + a a"^ + a 



(95) 



'2 sinh {huj /T) ^2 sinh {Tiuj/T) ^ 
with the density operator from eq. (^). For example, the mean number of original quanta 
is given by 

VV2 sinhr^y sinhr/ /iwsinhr/ Asinh^ihu /2T) ^ ' 

in agreement with ([771). The distribution of the original quanta differs from the Planck 
formula (0). At a narrow range of noise, T ^ huj, A ^ exp(— /;,u;/2T), which corresponds 
to the quantum low-temperature limit, we obtain a normal quantum result, 

k = e-^'^, T^huj. (97) 

A broad noise, A T /huj ^ 1, leads to the classical limit of high temperature. In this limit, 
the number of quanta is increasing with temperature quadratically rather than linearly, 

rp ^ 2 



T > huj. (98) 
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Let the random field applied to an oscillator with the electric charge e be abruptly 
removed at the moment t = 0. The oscillator remains excited and starts to radiate electro- 
magnatic waves losing its energy with the rate 

3c3 3c3 sinh(^u;/r) ~ Sc^ 4smh\hu/2T)' ^ ^ 

The emission rate directly measures effective temperature. Having no knowledge about the 
noise properties and detecting only the frequency u and the intensity E of radiation, we 
could assume thermal equilibrium inside the radiating system and assign temperature Tq 
by the direct application of the Planck formula instead of the actual distribution (|79D . It is 
easy to see that these temperatures are related through 

smh\hu/2T) = i smh{huj/2To)e^'^^^'^\ (100) 

For the narrow noise, two definitions agree. 

To ^ T, T < hiu. (101) 

However, under conditions of the broad noise, they lead to very different assignments with 
To>T, 

To ^ — , T > huj. (102) 

hio 

Two definitions can be discriminated in the case of the spectrum of normal modes with 
different frequencies by the deviations from the Planck formula at the low-frequency edge, 
if the noise amplitude F is not frequency dependent. 

VI. A MANY-BODY EXAMPLE 

Here we illustrate the concept of correlational entropy by calculating this quantity for 
a realistic many-fermion model with strong interparticle interaction of a variable random 
strength. We use for this purpose the version of the nuclear shell model which provides the 



best available description of spectroscopic data for all sd-imc\ei The model hamiltonian 
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consists of the one-body part given by a spherical mean field due to the closed shell core of 
^^O and the semiempirical residual two-body interaction determined by 63 reduced matrix 
elements {{jij2)JT\\V\\{j3j4)JT) allowed in this truncated Hilbert space by the conservation 
laws of angular momentum J and isospin T. Apart from numerous specific nuclear calcula- 
tions, this model was recently used p5| , pB| ,[7| for studying highly excited states beyond the 



limits of experimental information. One of the methods utilized in this analysis studies the 
evolution of observables as a function of the strength of the residual interaction [0,^,^]. 

We take as a generic example a set of states with quantum numbers J'^T = 0~^0 in the 
^^Mg nucleus with 4 protons and 4 neutrons in the srf-shell model space. The dimension 
of this set of states is =325. The unperturbed basis \k) is that of simple independent 
particle configurations which are projected onto correct values of integrals of motion J and 
T. The residual interaction is split into two parts, diagonal and off-diagonal with respect to 
the basis \k). The diagonal part is included into unperturbed energies in order to lift the 
degeneracy of bare configurations. The off-diagonal interaction with the overall factor A in 
front is our "random noise" with A = and A = 1 corresponding to the independent particle 
case and to the realistic strength, respectively. 

The many-body hamiltonian is diagonalized as a function of A which provides us with 
325 energy terms Ea{X). These terms were shown in Fig. 1 of Ref. |3^. Multiple avoided 



level crossings rapidly transform the eigenstates of the hamiltonian into complicated super- 
positions of many original configurations. The generic signatures of quantum chaos in level 
statistics, such as nearest level spacing distribution, spectral rigidity, and level curvature 
distribution, agree with the GOE predictions already at A ~ 0.3. The calculation of in- 
formation entropy in the original mean field basis shows that the evolution of complexity 
continues so that at A = 1 the states in the middle of the spectrum are close to the GOE 
limit of / = ln(0.48iV). 

Using the known eigenfunctions \a) of the system in the original basis \k), we find the 
density matrix p^fc/(A) of Eq. and construct the ensemble of such matrices for different 
values of A. The uniform averaging over A in the chaotic interval 0.3 < A < 1 gives 
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the density matrix p"^") of Eq. This matrix can be diagonahzed and its eigenvalues 

(occupation numbers) determine correlational entropy (P). 

Fig. 3 shows correlational entropy S*" of all 325 states ordered according to increasing 
energy. Although the fluctuations at this dimension are significant, nevertheless we see 
the same qualitative pattern as the one obtained with the help of information entropy 
23|,[7[] . The degree of mixing measured by correlational entropy evolves along the spectrum 



revealing the regular increase with excitation energy (in finite Hilbert space the statistical 
properties of eigenstates are symmetric with respect to the middle of the spectrum; in the 



thermal description |^ the center corresponds to infinite temperature with positive and 
negative temperature regions on the left and on the right, respectively). The behavior of 
information entropy and of single-particle entropy defined by the evolution of the single- 
particle occupation numbers is essentially the same. However, the absolute scales are very 
different. As seen in Fig. 3, the typical values of = exp[S'^"-'] for the most complicated 
states are close to 6 which is much lower than the localization length exp found with 
the use of information entropy. From the perturbative analysis of Sect. 3 we know that 
the number of effectively occupied eigenstates of the density matrix is related to the order 
of perturbation theory. Then can be interpreted in terms of a number of particle-hole 
excitons created by the spectral evolution. This point deserves to be studied more in detail. 



VII. CONCLUSION 

We have discussed some examples of quantum-mechanical systems under the conditions 
of stationary external noise which was defined by a distribution function of parameters 
controlling the interaction hamiltonian. The maximum of quantum information on such 
systems is provided by the density matrix which was found for various cases, analytically in 
Sect. 3-5 and numerically in Sect. 6. The perturbative treatment of Sect. 3 illustrates the 
mechanism of equlibration through the consecutive selection of directions in Hilbert space 
which form the eigenbasis of the density matrix. In other examples the weakness of the 
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interaction was not assumed, and the solution was exact. 

We have used correlational (von Neumann) entropy as a tool for measuring the degree of 
disorder and complexity of stationary states. By definition, this entropy vanishes for pure 
quantum states. Therefore it characterizes the system in a given noise ensemble when the 
stationary states are mixed. In general, the arising configuration with certain occupation 
numbers of the eigenstates of the density operator does not coincide with thermal equilibrium 
as defined by canonical Gibbs ensembles. Accordingly, correlational entropy is not equal to 
thermodynamic entropy. Moreover, correlational entropy is calculated for the individual 
energy terms which evolve adiabatically as a function of the noise strength. The resulting 
state is in general different for different terms. 

It would be interesting to determine the conditions for the external noise which would 
give the same equilibrated state as in the heat bath. It is assumed usually that the necessary 
ingredient is the continuous spectrum of the normal modes represented in the spectral ex- 



pansion of noise . We deal with the stationary noise represented by random parameters 
in the hamiltonian. In Sect. 5 we have however shown that the ground state of a harmonic 
oscillator in a random uniform field can be described as an equilibrated thermal state of an 
effective oscillator with temperature determined by the Gaussian width of the field distri- 
bution. The situation here is similar to that in many-body physics where the interacting 
system can be modeled |Q by a gas of dressed quasiparticles with properties depending 
on temperature (in our case the energy spectrum of the oscillator is not renormalized but 
the coordinate scale is determined by noise and changes with effective temperature). The 
difference between a simply heated oscillator in a thermostat and an oscillator excited by a 
Gaussian noise and described with the aid of effective temperature might be important for 
the problems of multiple meson production P6| , |5T| in high energy collisions. We can remind 



also that the notion of the effective quantum oscillator appears naturally in the problem of a 



uniformly accelerated observer in the Minkowski world ||5^J5^ . An observer falling with the 
proper acceleration g sees the Minkowski vacuum as a black body emitter with the effective 
temperature T = hg/{2TTc). This consideration is closely related to the Hawking black hole 
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radiation ||5^. However, in those cases, one has squeezed rather than coherent states of the 
oscillator (in terms of the unperturbed system they are produced by the source creating the 
quanta pairwise which is described with the help of the Bogoliubov transformation). 

Representing the complexity of individual quantum terms, our invariant correlational 
entropy can be juxtaposed to representation-dependent information entropy. We drew the 
attention of the reader to their similarity and distinction in various applications. Although 
they are formally quite different and may react in a different way on the change of parame- 
ters of a simple regular system. Sect. 4, they behave qualitatively similar in a complicated 
system. Sect. 6, where one sees the standard signatures of quantum chaos. Due to similarity 
of adjacent states in the chaotic regime 0, both entropies are smooth functions of excita- 
tion energy and therefore can be considered as thermodynamic variables. Both entropies 
carry information on complexity of individual states and its evolution along the spectrum. 
Quantitatively, information entropy (in the appropriate basis!) expresses this complexity in 
terms of a number of mixed simple configurations whereas correlational entropy measures 
essentially similar properties in larger blocks as a number of classes of states effectively 
mixed by external noise. Of course, one should remember that information entropy refers 
to a given hamiltonian while correlational entropy describes a "system plus noise" complex. 
The further studies will bring the more deep insight into the problem. 
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Figure captions 

Figure 1. Correlational entropy of the two-level system (16) as a function of the ratio 
T = IV I k of the strength of the mixing interaction 2V to the range of the random fluctua- 
tion A of the level positions; parts a,b,c and d correspond to the values x = e/A = 0.2, 0.6, 1 
and 5, respectively, where e is the static level spacing. 

Figure 2. Correlational entropy of the two-level system (16) as a function of the ran- 
dom perturbation strength for four values of the intrinsic mixing strength C, = t jx'- 
0.01 (sohd hue), 0.1 (dots), 2 (dash), 10 (dot-dash). All the curves asymptotically approach 
the limiting value In 2; the rate of approach is slower for larger values of the strength pa- 
rameter. 

Figure 3. Correlational entropy of 325 J'^T — 0+0 states in the sd-shell model for the 
^''Mg nucleus. The states are ordered in increasing energy. 
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